Bragg spectroscopy of the multi-branch Bogohubov spectrum 
of elongated Bose-Einstein condensates 
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We measure the response of an elongated Bose-Einstein condensate to a two-photon Bragg pulse. 
If the duration of the pulse is long, the total momentum transferred to the condensate exhibits 
a nontrivial behavior which reflects the structure of the underlying Bogoliubov spectrum. It is 
thus possible to perform a spectroscopic analysis in which axial phonons with different number of 
radial nodes are resolved. The local density approximation is shown to fail in this regime, while 
the observed data agrees well with the results of simulations based on the numerical solution of the 
Gross-Pitaevskii equation. 



PACS numbers: 

The crossover between phonon and single-particle ex- 
citations in the Bogoliubov spectrum of a weakly inter- 
acting Bose gas jlj is one of the main "textbook" con- 
cepts that can be directly tested in the case of trapped 
Bose-Einstein condensed gases. Bogoliubov quasiparti- 
cles have been already produced in elongated conden- 
sates by using two-photon Bragg scattering |H, ||, ^ ||. 
For excitations with frequency lo and wavevector k, the 
phononic character has been checked by observing that 
the static structure factor S{k) is less than 1 [||, by 
measuring the quasiparticle amplitudes Uk and W/t ||^, ^ , 
and by showing that the dispersion relation w(fc) is linear 
at low fc [g |. 

In all of these cases, the local density approximation 
(LDA) has been used to adapt the Bogoliubov theory of 
uniform gases to the actual inhomogeneous condensates. 
This approach is expected to be accurate for large con- 
densates, where the density profile varies smoothly on 
the scale of the excitation wavelength and the system 
behaves locally as a piece of uniform gas with a local 
Bogoliubov spectrum |^. This applies, for instance, 
to axial excitations of elongated condensates, when the 
wavelength of the excited states is much smaller than the 
size of the system along the major axis. These excita- 
tions can then be classified with a continuous wavevector 
k. However, the finite transverse size of the condensate 
also produces a discreteness of the spectrum, which is 
ignored in LDA (see |^ for a recent classification of nor- 
mal modes in anisotropic condensates and [|lO[ |l^ for the 
limiting case of an infinite cylinder). 

In this Letter we show that the response of the con- 
densate to a Bragg pulse is indeed significantly affected 
by the radial degrees of freedom. In particular, if the 
duration of the pulse is longer than the radial trapping 
period, the condensate responds resonantly at the fre- 



quencies ujn^{k) of axial quasiparticles with nodes in 
the radial direction. By using Bragg pulses longer than 
in previous measurements we resolve this multi- 

branch spectrum, finding good agreement with the pre- 
dictions of Gross-Pitaevskii (GP) theory. 

As described in Ref. our condensate consists of 
N — 10^ atoms of ^^Rb, with a thermal fraction of 5% 
or less. The radial and axial trapping frequencies are 
uj±_ = 27r(220Hz) and uj^ = 27r(25Hz) respectively. The 
radial and axial Thomas-Fermi radii of the condensate 
are R = 3.1^m and Z = 27.1/im. 

We excite the condensate by using two Bragg beams 
with approximately parallel polarization, separated by an 
angle 9. The Bragg beams illuminate the entire conden- 
sate for a time Ib- The beams have a frequency differ- 
ence LO determined by two acousto-optic modulators. If a 
photon is absorbed from the higher-frequency beam and 
emitted into the other, an excitation is produced with 
energy hio and momentum hk, where k = 2kpShi(0/2), 
and kp is the photon wave number. The wave vector k 
is adjusted to be along the z-axis. To measure a single 
point on the excitation spectrum cj(fc), k is fixed by 9, 
and uj is varied. 

The measured quantity is the total momentum trans- 
ferred to the condensate along z. This is obtained by 
switching off the trapping potential after the Bragg pulse 
and taking absorption images of the density distribution 
of the expanding condensate. In Fig. 1 we show typical 
results for short Bragg pulses of duration = 1 msec. 
The observed momentum is plotted as a function of 
UJ for three different values of fc. In Ref. [Q we made 
Gaussian fits to these type of curves, taking the position 
of the maximum to be the quasiparticle dispersion uj{k) 
and showing that this dispersion was consistent with the 
Bogoliubov spectrum in the LDA. 
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FIG. 1: Pz, in arbitrary units, as a function of uj for short 
pulses (ts = 1 ms). Bottom, center, and top plots correspond 
to different fc-values. Points with error bars are the measured 
values. Solid lines are the results of GP simulations. The 
dashed lines are the LDA predictions. 

A deeper analysis of the same results can be performed 
by numerically solving the GP equation for the order 
parameter of the condensate $(r, i) 

indt^=(-^ + V + gm^U (1) 

where g = Airh^a/m is determined by the s-wave scatter- 
ing length a. The external potential can be taken as the 
sum of the harmonic confinement and the Bragg poten- 
tial, 

Tfi 

V{T,t) = —{ujlr\+ulz'^) + e{t)VBCOs{kz^Lot) (2) 

where r^—x^^ and 6{t) is equal to 1 in the interval 
< t < ts and outside. 

The ground state at i = can be found as the station- 
ary solution of Eq. (|l]). Then, the time dependent GP 
equation can be solved at t > to simulate the Bragg 
process |l^ . We take advantage of the axial symmetry 
of V to map the order parameter into a two-dimensional 
grid of points N± x (typically, 64 x 1024) and evolve it 
by means of a Crank-Nicholson differencing method with 
alternating direction implicit algorithm, as in [ p^ . 

The momentum transferred to the condensate can 
be calculated from $ through the definition — 
{—ih/2) J dr $*9z$ -I- c.c. Fig. 1 shows good agreement 
between the results of the GP simulations (solid lines) 
and the experimental data, for short pulses. The strength 
Vb in the GP equation is used as a free parameter, so the 



comparison is restricted to the position and shape of the 
peak. In principle, Vb could be determined by using 
the experimental parameters for the power, area, and di- 
rection of polarization of the two laser beams, but this 
estimate might be affected by a significant uncertainty. 

The LDA curves (dashed lines in Fig. I), which are 
obtained by using the local Bogoliubov spectrum to cal- 
culate Pz{k,Lu) in the linear response regime as in are 
also reasonably close to the GP predictions. For shorter 

we checked that the accuracy of the LDA is even bet- 
ter, as discussed in ^. 

A rather different situation is found for longer pulses, 
with of the order of the radial trapping period (about 
4.5 msec) jl^. Fig. 2 shows the momentum transferred in 
the GP simulations (solid curves) for various time dura- 
tions tB- The lowest curve, after 1 msec, is a broad peak 
as in Fig. I, close to the LDA (dashed curves). For later 
times however, a multi-peak structure appears, strongly 
deviating from LDA. We find a similar behavior for all 
values of k in our simulations, and for a wide range of 
intensities Vb- 
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FIG. 2: The time dependence of the momentum transferred 
during the Bragg pulse, for k — 3.f/^m~^. Solid curves are 
the results of GP simulations with Vb = Q.2hu}±_ plotted at 
various times. Starting with the lowest solid curve, Ib = 
1, 2, 3, 4, 5 and 6 msec. Dashed curves are the corresponding 
LDA predictions. The arrows indicate the GP predictions for 
the frequencies of the normal modes of the condensate, having 
0, 1 and 2 radial nodes and the same axial wave vector k. 

Figures 3a and 3b show measurements of for k = 
1.4/im~^ and 3.1/xm~^ respectively, with long pulses of 
duration 10 msec and 6 msec, respectively. The intensi- 
ties of the Bragg beams are adjusted so that the num- 
ber of excitations created on resonance is no more than 
roughly 25% of the number of atoms in the condensate. 
Each point in Fig. 3 is an average of about 5 measure- 
ments. The results of GP simulations are indicated by 
the upper solid curves, which are obtained by choosing 
Vb = 0.2?ki;_L and 0.35Sw_l for Figs. 3a and 3b, respec- 
tively. Both the measurements and the GP simulations 
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show multiple peaks, due to the excitations of quasipar- 
ticles with different n^. 

In Fig. 3, we see good agreement between the locations 
of the peaks in the experimental and GP results, but the 
experimental peaks are broader due to noise. The ma- 
jor source of noise is the sloshing of the condensate in 
the trap. Specifically, sloshing of speed Vsi in the axial 
direction yields a Doppler shift in the light potential of 
frequency kvsi/^iT. This gives a frequency noise that sig- 
nificantly broadens the peaks of P^. However, we can 
reduce this noise by estimating the velocity of the con- 
densate in the trap at the end of the Bragg pulse from 
its position in the time-of-flight image, and by explicitly 
adding the corresponding Doppler shift to the applied u 
in the laboratory frame. This correction requires that the 
Bragg pulse be much shorter than the axial trap period 
(40 msec). This requirement is sufficiently met by the 6 
msec pulses employed for k = 3.1/xto~^. Therefore, the 
data points of Fig. 3b contain the correction. 

We verify that the multi-peak structure reflects the 
fundamental normal modes composing the multi-branch 
spectrum, and does not depend significantly on the inten- 
sity Vb- The lower curves of Figs. 3a and 3b correspond 
to the same simulations as the upper curves, but with 
Vb reduced by a factor 10. In this case, the system is 
excited in the linear response regime. The momentum 
transferred is two order of magnitude smaller, but the 
resulting curves have peaks at the same locations as for 
larger Vb- Furthermore in the case of A: = 3.1/im~^ we 
repeat the measurement and simulation of Fig. 3b with 
an intensity 30% greater and find peaks again at the same 
location. This means that nonlinear effects are not cru- 
cial in our observations. Nevertheless, they might be 
interesting. Indeed, looking at Fig. 3a one notices small 
differences between the lower and upper curves, for small 
and large Vb respectively. The curve for large Vb dis- 
plays two side peaks around the main peak at 0.4 kHz. 
These side peaks, whose shape significantly depends on 
ts, are not visible in the linear regime of the small Vb 
curve and might be due to nonlinear effects. Similar ef- 
fects were found by a one-dimensional simulation in Ref. 

Now we use the GP equation to directly find the multi- 
branch Bogoliubov spectrum of the condensate. Specifi- 
cally, we analyze the oscillations in the condensate den- 
sity induced by the Bragg process. In the GP simulation 
we let the condensate freely oscillate in the trap after the 
Bragg pulse and we perform a Fourier analysis of the den- 
sity variations. This analysis shows that for each k and w, 
the density oscillates as a superposition of modes of fre- 
quency ujn^{k), which are excited by the Bragg potential 
in Eq. (^) due to the inhomogeneity of the condensate in 
the radial direction. For symmetry reasons, only modes 
with azimuthal angular momentum m = are excited. 
The calculated frequencies are shown as open circles in 
Fig. 4 ijl^ . The lowest branch corresponds to Bogoliubov 
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FIG. 3: Pz in units of Nhk, as a function of lj for long Bragg 
pulses. For (a) and (b), t_B = 10 msec and 6 msec respectively, 
and k — 1.4fim~^ and 3.1^m~^ respectively. The circles are 
the measured values. The upper curve is the result of GP 
simulations with Vb = O.S5huj± and 0.2huj±, for (a) and (b) 
respectively. The lower curve is the same but with Vb ten 
times smaller. For the lower curve, Pz is multiplied by 10 
and 20 for (a) and (b) respectively. The arrows on the uj axes 
indicate the normal mode frequencies, as in Fig. 2. 

axial modes with no radial nodes. The second branch has 
one radial node; it starts at 2lu± for k = 0, where it cor- 
responds to a purely radial breathing mode. The third 
branch has two nodes |lj] . 

In the limit of small oscillations, all these states co- 
incide with the solutions of the linearized GP equation, 
i.e., the generalization of the Bogoliubov equations to in- 
homogeneous condensates. For an infinite cylinder these 
solutions were calculated in Ref. jlOj and, in the hydro- 
dynamic limit (fc much smaller than the inverse healing 
length ^^^), in Ref. |]ll]. For k <^ our multi-branch 
spectrum turns out to be very close to the spectrum pre- 
dicted in Ref. Q, indicated by dashed hues in Fig. 4 
(also see Fig. 1 of Ref. [0). is indicated by a dotted 
line in Fig. 4. 
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the Gross-Pitaevskii equation. The local density approx- 
imation fails to reproduce these multiple radial modes. 

The technique of long Bragg pulses could also be used 
to resolve inherent dissipation mechanisms, such as in- 
homogeneous broadening and coupling between modes, 
including Beliaev and Landau damping. The tech- 
nique could be particularly effective if the inhomogeneous 
broadening could be overcome, by means of echo tech- 
niques, or the use of non-quadratic optical dipole traps. 
Finally, very sensitive GP simulations could be used in 
combination with long pulses, to carefully study non- 
linear dynamics. 

This work is partially supported by MIUR-COFIN2000 
and the Israel Science Foundation. 
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FIG. 4: Excitation frequencies vs. k for a trapped BEG. The 
open circles are the frequencies of the normal modes with rir — 
0, 1, 2, . . . nodes in the radial direction, as calculated from the 
GP equation. The filled circles are the position of the maxima 
of the main peaks in Pz for the long-pulse measurements of 
Fig. 3. The dashed lines are the prediction for an infinite 
cylinder with k <C C^^. The dotted line indicates the inverse 
healing length. 



The calculated normal mode frequencies are shown as 
arrows on the oj axes in Figs. 2 and 3. Both the theoreti- 
cal curves and the experimental results clearly show that 
the momentum is resonantly transferred to the conden- 
sate when the Bragg frequency is close to the Bogoliubov 
branches. 

The locations of the measured peaks of Fig. 3 are also 
shown as filled circles in Fig. 4. The good agreement be- 
tween these measured peaks and the simulated spectrum 
(open circles) identifies the peaks in the observed spec- 
trum as being axial quasiparticles with n,- radial nodes. 
Bragg pulses can indeed be used to perform high resolu- 
tion spectroscopy of the multi-branch Bogoliubov spec- 
trum of trapped condensates in a regime where the LDA 
is not applicable. 

In conclusion, we have used long Bragg pulses to spec- 
troscopically measure the first and zeroth order radial 
modes in the phonon spectrum of a BEG. These high- 
resolution measurements agree well with simulations of 
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A longer timescale is fixed by the axial trapping pe- 
riod (40 ms). For such durations other processes be- 
come important (reflection of axial modes at the bound- 
ary, center-of-mass motion, nonlinear effects, etc.). This 
regime is considered in WA. 
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The calculated frequencies have an uncertainty of the or- 
der of 0.05 kHz, due to the finite time interval used for 
the Fourier analysis. 

In this range of k only three branches are visible in the 
Fourier spectrum. The efficiency of the Bragg perturba- 
tion (Q) in exciting higher branches increases with k. 



